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Abstract 

We propose a method to obtain an improved Hamihonian (action) for 
the Ising universahty class in three dimensions. The improved Hamil- 
tonian has suppressed leading corrections to scaling. It is obtained by 
tuning models with two coupling constants. We studied three different 
models: the ±1 Ising model with nearest neighbour and body diagonal 
interaction, the spin-1 model with states 0, ±1, and nearest neighbour 
interaction, and (/>^-theory on the lattice (Landau-Ginzburg Hamilto- 
nian). The remarkable finite size scaling properties of the suitably 
tuned spin~l model are compared in detail with those of the standard 
Ising model. Great care is taken to estimate the systematic errors 
from residual corrections to scaling. Our best estimates for the criti- 
cal exponents are v = 0.6298(5) and r/ = 0.0366(8), where the given 
error estimates take into account the statistical and systematic uncer- 
tainties. 
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1 Introduction 



In Monte Carlo simulations the size of the systems that can be studied is lim- 
ited by the memory of the computer and by the CPU-time that is available. 
Therefore, in many instances, finite size scaling is the key to a precise 
determination of properties of statistical systems at criticality. Finite size 
scaling laws are affected by corrections to scaling. These corrections to scal- 
ing cause systematic errors in the results for the infinite volume limit one is 
interested in. With improving statistical accuracy of the Monte Carlo data it 
becomes important to deal properly with systematic errors. One way to pro- 
ceed is to include corrections to scaling into the fit ansatze when analysing 
the data. Another, more fundamental way is to remove corrections already 
from the system to be studied. 

Renormalization group (RG) [0] offers (at least in principle) a way to 
achieve this goal. RG fixed point actions^ are free of corrections to scaling. 
However such actions in general contain an infinite number of couplings. In 
practical applications one is forced to truncate the action to a finite number 
of terms, which in fact is an uncontrolled approximation. For an applica- 
tion of this strategy to asymptotically free models see the work on Perfect 
Actions [§. 

A different approach was pioneered by K. Symanzik [Q. Higher order 
terms are added to the action. By imposing certain conditions on observables 
leading corrections to scaling are eliminated. While Symanzik formulated his 
method in the framework of perturbation theory, recently there have been 
attempts to apply this method in a non-perturbative, i.e. numerical, setting 

i- 

Our present approach is closer to this latter point of view than to the 
block spin renormalization group inspired framework. 

The idea to improve the scaling properties of Ising models by moving 
to models with generalized actions and tuning the coupling constants as to 
obtain reduced corrections to scaling was already followed in the work of 
Blote et al. 0. However, in their work they do not uncover the principle and 
method they used to obtain their improved actions. 

In a recent paper 0, we presented data obtained with a spin-1 Ising 
model with states 0, ±1. Tuning the two coupling constants in the proper 

*What is called Hamiltonian or energy function in Statistical Mechanics is called Action 
in Euclidean Field Theory. 



4 



way, we were able to reduce the corrections to scaling in various quantities 
dramatically. Especially the Binder cumulant and its derivatives, and also the 
susceptibility, could be fitted with scaling laws without corrections to scaling 
terms, yielding very precise estimates of the 3D Ising critical exponents. 

Recently the authors of refs. P, |^ argued that the improvement discussed 
above does not lead to reduced error estimates for critical exponents. They 
argue that the error estimates given in our recent paper [0] are underes- 
timated, since we do not take into account residual leading correction to 
scaling corrections. Such corrections might well be present, since the param- 
eters of our improved model are computed numerically. However, this is not 
the full story as we shall explain in this paper. Our argument is based on 
the fact that ratios of correction to scaling amplitudes are universal. 

In this article, we describe in detail our method, the numerical results, 
and the fitting procedures. We confront the results from the improved ac- 
tions with high precision data from simulations of the standard Ising model, 
estimating with a well defined procedure systematic errors for both actions. 

2 Improving the Scaling Behaviour 
2.1 The Models 

Usually, Monte Carlo studies of the Ising model are done using what in field 
theory is called 

Standard Action 

S = -f3j2s.s,. (1) 

<i,j> 

The Si take values ±1, and the spin-spin interaction is a sum over all nearest 
neighbour pairs < i,j >. A precise estimate for the critical coupling was 



obtained in reference 0: /^^ = 0.2216544(3) (3). 

In the following we will introduce and study three different models in the 
3D Ising universality class, each of them governed by two coupling constants. 
In all three cases the Boltzmann factor is given by exp(— S"). 

Spin-1 Model 

S = -/3 Y: s,s, + DY.sI (2) 

<i,j> i 
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The Sj take values 0, ±1, and the spin-spin interaction is a sum over all nearest 
neighbour pairs. This model was introduced and studied in 0. There, D 
was fixed to In 2. The critical /3 corresponding to this particular value of D 
was estimated in § to be 13^ = 0.3934224(10). 

NNN Model 

S = -pi J2 '^i '^j - P2^Si Sj . (3) 

<i>j> 

The Sj take values ±1, and the spin-spin interaction is a sum over all nearest 
neighbour pairs < i,j > and third neighbour pairs (body diagonals) 
B15te et al. fixed P2/P1 = 0.4 and obtained /?i,c = 0.1280036(5). 

04-Model 

s=-(3y: ^■^'P^ + e + a e(</'' - 1)' • (4) 

<i,j> i i 

The variables assume real values. In the limit A ^ 00 one recovers the 
standard Ising model. 

The three models defined above have in common that they have a second 
order critical line in the space spanned by the two coupling constants. We 
exploit the degree of freedom of moving on the critical line to find critical 
models with reduced corrections to scaling. How this is done will be explained 
in the following subsections. 



2.2 Matching of Phenomenological Couplings 

We study two independent phenomenological couplings of the 3D Ising model, 
to be called i?j, z = 1,2 in the following. Both quantities are universal, i.e. 
at criticality their infinite volume limit R* does not depend on details of the 
microscopic Hamiltonian. Ri is the ratio of partition functions with anti- 
periodic and periodic boundary conditions, respectively, 

Rl = Za/Zp. (5) 

The lattices will always be cubical with extension L in each of the three 
directions. Antiperiodic boundary conditions are imposed only in one of the 
three lattice directions. R2 is the Binder cumulant, 

^2 = Q = ^. (6) 
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Here, m denotes the magnetization per spin, 

i 

The Ri are L-dependent and, of course, functions of the couphng parameters 
in the action. For the two-couphng models defined above, we define "flows" 
(hues of constant physics) {Ki{L), K2{L)) by requiring that 

R,mK,{L),K2{L)) = Rl. (8) 

Ki and K2 represent the two coupling constants of the model. In the next 
subsection wc shall demonstrate that with increasing L the flows of (A'l, K2) 
converge towards a critical point which has no leading order corrections to 
scahng. 



2.3 RG Analysis of the Matching Condition 

The main features of the two-couphng models can be discussed in the frame- 
work of the renormalization group. The scaling properties can be derived 
from the linearized RG transformation at the fixed point. 

We consider general Hamiltonians with couplings K^., where a = 1, 2, . . . 
An RG transformation, realized, e.g., by a block spin transformation, changes 
these couplings according to 

K ^K' ^ R{K) = K'{K) . (9) 

A fixed point K* is defined through R{K*) ~ K*. The linearized transfor- 
mation at the fixed point can be represented by a matrix 



(10) 



K=K* 



One introduces "normal coordinates" (scaling fields) by 

u,^u,{K)^Y.^iAKa-K), (11) 

a 

where </7j denotes the i-th (left) eigenvector of the matrix T, 



7 



The Ui transform under RG transformations like 

Ui XiUi . 

In Ising type models the leading eigenvalues are given by 

\i = h"\ A2 = 6-^ \-i = b-\ (12) 

where x » uo. Note that ,x = 2 for the lattice Gaussian model. From leading 
order e-expansion one expects that x is close to 2 at the Wilson-Fisher fixed 
point, h denotes the scale factor of the RG transformation. 

Let us now assume that we have only two non-vanishing couplings Ki 
and K2 in our Hamiltonian. Let us then write down explicitly the condition 
for being critical {ui = 0) and eliminating the leading corrections to scaling 
(^2 = 0). The first condition reads 

- Kl) + (^i,2(X2 - = Ki,3 , (13) 

whereas the condition U2 = translates to 

(^2,1(^1 - K*) + ^2,2(^2 - K;) = K2,s . (14) 
For i = 1,2, the Kj^a are given by 

a>3 

Let us now study how our matching procedure with the two quantities 
-Ri = Za/ Zp and R2 = Q works. The are functions of the bare couplings 
and the lattice size: 

Rk = Rk{L,K,,K2). (16) 
We express these quantities as functions of the scaling fields defined above, 

Rk{L, Xi, K2) = Rk L--^'^) . (17) 

Here, the upper index indicates that the scaling field is taken at the scale 
of the lattice spacing. The pref actor promotes the scaling field to its value 
at the scale L. Taylor-expansion of the Rk around their fixed point values 
yields 

Rk^Rl + Tk^L^'^'u^l^ + Tk,2L-^U^2^ . (18) 
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The matching conditions Rk = Rl are thus equivalent to 

rfc,iLi/'^4')+rfc,2L--4') = 0, (19) 

for k = 1,2. We obtain, as desired, the solution u^i^ = (criticality) and 
'^2^'' = (iio leading order corrections). Including higher order corrections in 
the scaling ansatz, governed by and exponent A3, one can convince oneself 
that fixing Ri and R2 to their fixed point values leads to convergence to the 
critical line Ui = with corrections that decay like L^^"^/" . The M2 = 
condition is approached with a much slower rate, namely like L'^^"^ . 

2.4 Computing the Matching Flows 

For the three 2-coupling models specified above, we set up a procedure to de- 
termine the flows of couplings {Ki{L), K2{L)) such that eq. (H) was fulfilled. 
To this end we used estimates for Rl = 0.5425(10) and R*2 = 0.6240(10), 
which were an outcome of a preliminary analysis of data obtained using the 
standard action. 

The matching couplings were searched for using a Newton iteration, based 
on the inversion of a matrix made up from Monte Carlo estimates of the 
derivatives of the Ri with respect to the two couplings. Typically three 
to four iterations were sufficient to find couplings such that Za/Zp and Q 
attained the prescribed values within the given statistical precision. The 
results are given in table I. 

A first look at the table reveals that both for the spin-1 and the 0^-model, 
the flow converges to a fixed point quickly, whereas it keeps moving strongly 
in the case of the NNN model. A rough explanation of this could be the 
following: In order to remove the leading corrections to scaling, one has to 
move efficiently between the Gaussian model (equivalent to the 0^-model with 
vanishing A) and the non-trivial Wilson-Fisher fixed point. Moving between 
these two fixed points is most efficiently done using a 0^-type coupling, which 
is implicitly present also in the spin-1 model. The NNN model seems to 
need renormalization to a larger scale in order to come close to the flow line 
connecting the Gaussian with the Wilson-Fisher fixed point. 

Plotting the second coupling vs. the first one, one finds that with very 
good precision the critical line can be approximated by a straight line: 

K2{L) = ai + a2 Ki{L) , (20) 

with 
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Spin-1 Model 



L 


/3 


D 




Q 


3 


0.35737 


0.4401 


0.54201(30) 


0.62447(21) 


4 


0.37250 


0.5510 


0.54242(25) 


0.62347(18) 


5 


0.37794 


0.5883 


0.54292(19) 


0.62421(15) 


6 


0.38210 


0.6169 


0.54221(49) 


0.62426(35) 


7 


0.38419 


0.6311 


0.54326(30) 


0.62366(22) 


8 


0.38320 


0.6241 


0.54291(37) 


0.62430(27) 


9 


0.38320 


0.6241 


0.54259(60) 


0.62403(44) 


10 


0.38320 


0.6241 


0.54288(57) 


0.62431(42) 






NNN Model 




L 






Za/ Zp 


Q 


4 


0.12266 


0.05406 


0.5429(2) 


0.6238(2) 


5 


0.12928 


0.05028 


0.5427(1) 


0.6241(2) 


6 


0.13431 


0.04734 


0.5425(1) 


0.62441(8) 


7 


0.13800 


0.04518 


0.5425(3) 


0.6242(2) 


8 


0.14069 


0.04361 


0.5426(2) 


0.6243(1) 


9 


0.14292 


0.04231 


0.5430(2) 


0.6232(1) 


10 


0.14406 


0.04165 


0.5425(2) 


0.6242(2) 


11 


0.14590 


0.04059 


0.5429(1) 


0.62385(8) 


iz 


0.i47z4 


0.03981 


0.5426(1) 


0.6z35(ij 


13 


0.14808 


0.03933 


0.5431(1) 


0.6236(1) 








Model 




L 




A 


Za/Zp 


Q 


3 


0.35303 


1.5248 


0.5420(3) 


0.6244(2) 


4 


0.36338 


1.3282 


0.5425(3) 


0.6243(2) 


5 


0.36908 


1.2188 


0.5425(2) 


0.6241(1) 


6 


0.37165 


1.1689 


0.5427(2) 


0.6241(1) 


7 


0.37270 


1.1481 


0.5424(2) 


0.6243(1) 


8 


0.37308 


1.1410 


0.5421(1) 


0.6245(1) 


9 


0.37273 


1.1479 


0.5416(3) 


0.6247(2) 



Table 1: Flows of couplings defined such that for all lattice sizes 
the two quantities Ri = Za/Zp and R2 = Q match (to the given 
statistical precision) with their fixed point values Rl = 0.5425(10) 
and R^ = 0.6240(10). 
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model oi 02 



-2.04 
0.1253 
8.29 



spin-1 
NNN 

0^ 



6.95 
-0.5804 
-19.17 



It is assuring that for the spin-1 and the NNN model, the critical coupling 
estimates determined by Blote et al. are in good agreement with our critical 
lines. 

Fits of Ki{L) with a power law 

iri(L) =ci + C2L-" (21) 

yielded good fits in all cases, with exponents of order 3 in the cases of spin-1 
and 0^-models. Note that this exponent is much larger than —x + uj that is 
to be expected from the theory presented above. For the NNN model, the 
exponent is of order 1. 

The slow convergence of the NNN flow motivated our decision to discard 
this model from further investigation. Also, the behaviour of the 0^ and 
spin-1 models appears to be very similar. We decided to concentrate on the 
spin-1 model, and leave the 0^-model for later study. 

Note that our result for the optimal A ~ 1.145 of the 0'^-model is consis- 
tent with the observation of ref. |]^ that the optimal A should be close to 
one. 



2.5 Identification of the W2 = Line 

After the decision to concentrate on the spin-1 model, we determined an 
approximation of the M2 = manifold. This was done by looking at the 
derivatives of the Ri at criticality. From the RG analysis of section p.3| one 
infers 

dRk ^ ORk dui , , _ 

The left hand side of the equation can be determined by Monte Carlo sim- 
ulation and will therefore be assumed as known. Taking into account that 
without loss of generality one can set 1 = 1, the equations can be solved 
(or fitted) if the left hand side is known at least for two different lattice sizes. 
We performed Monte Carlo simulations at (/5, D) = (0.3832, 0.6241) for lat- 
tice size 8, 9, 10, and 11. Fixing the exponents u = 0.63 and u = 0.81 [|lll we 
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obtained that the scahng field 02,2 should be approximately —1/3. Plugging 
this into the eq. (0), 

ip2^{K^ - Kl) + ¥^2,2(^^2 - = K2,3 , (23) 

one obtains that 3(3 — D should be kept constant. We used our simulation 
point (0.3832,0.6241) to fix this constant. In conclusion, one should approach 
criticality by varying /5 while adjusting D according to 

D = 3{p- 0.3832) + 0.6241 . (24) 

Of course, a precise estimate of Pc along this line still has to be determined. 
This will be discussed in the next section. 

The fit result for 2 can be used to perform a consistency check. In 
terms of the scaling fields, the critical line is characterized by 

iP - (3*) + ¥>i,2{D - D*) = , (25) 

ignoring all higher couplings. Solving for D and plugging in the "experimen- 
tal values" (3* = 0.3832, D* = 0.6241 and ipi,2 = -0.1439 one obtains with 
good precision the critical line approximation D = —2.04 + 6.95/3. Let us 
close this section by remarking that errors in the precise estimation of the 
^2 = line affect results for the critical exponents only weakly, e.g., the effect 
on the exponent u is of order 2,-1/'^-'^. 



3 Simulation Parameters and Statistics 
3.1 Standard Action 

Monte Carlo simulations of the 3D Ising model with standard action were 
performed at (3 = 0.2216545, which is a good approximation of the critical 
coupling |10 . 



For cubical lattices of size L = 2 up to L = 19 we performed simulations 



with the multi-spin demon update [^2|, |T3[. The update algorithm is local. 
Hence one expects a dynamical critical exponent z ^ 2. However, due to 
the multi-spin coding implementation a single sweep can be done substan- 
tially faster than with the cluster algorithm. Therefore one expects a better 
performance for the demon update on such small lattices. 
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L 


M 


S/V 


TE 




n 


6 


300000 


0.38412(29) 


1.12(2) 


1.18(2) 


0.53(1) 


8 


300000 


0.33222(27) 


1.13(2) 


1.19(2) 


0.56(1) 


12 


300000 


0.27005(24) 


1.13(2) 


1.18(2) 


0.57(1) 


16 


300000 


0.23277(23) 


1.17(2) 


1.19(2) 


0.58(1) 


24 


300000 


0.18866(20) 


1.15(2) 


1.13(2) 


0.55(1) 


32 


700000 


0.16214(13) 


1.18(2) 


1.14(2) 


0.56(1) 


48 


450000 


0.13098(14) 


1.20(2) 


1.09(2) 


0.54(1) 


64 


600000 


0.11258(11) 


1.26(2) 


1.10(2) 


0.55(1) 


96 


340000 


0.09075(13) 


1.26(4) 


1.05(3) 


0.53(1) 



Table 2: Testing the performance of the wall cluster algorithm. 
M denotes the number of measurements, separated by a single 
wall cluster step. The r's are integrated autocorrelation times. 
The unit of time is set to one complete update of the lattice; i.e. 
one measurement /(S*/!^). S/V denotes the average sum of the 
sizes of the flipped clusters, normalized to the lattice volume L^. 



For a subset of the smaller lattices, and for bigger lattices up to size 128, 
cluster update was performed, using a new variant of the algorithm, the wall 
cluster algorithm. Note that there is quite a lot of freedom in the selection 
of clusters which are flipped during one update step. In the wall cluster 
algorithm, one flips all clusters that intersect with a given lattice plane. 
Sequentially one takes lattice planes in 1 — 2 , 1 — 3 and 2 — 3 direction. The 
position is chosen randomly. Let us call the procedure to generate and flip all 
the clusters connected to the selected plane a wall cluster update step. The 
motivation for choosing this type of update was that the construction of all 
clusters that have elements in a lattice plane is needed for the measurement 
of Za/Zp anyway. Testing the performance of the algorithm shows that there 
is some small gain in efficiency compared to the single cluster update |T^ . 



Results for cluster size distribution and relaxation times of simulations at (3c 
are given in table ^. 

The average sum of the sizes of the clusters per volume that are flipped 
in one step behaves as 

S/V = CL'' . (26) 

Fitting the data of table |^ to this law, discarding the results from L < 24 
yields S/V = 1.008(4)L-°-527(i) ^^2^^^^^Qj^y rj.^^ integrated autocorrela- 
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tion times of the energy, the susceptibihty, and of the Za/Zp measurements 
(see below) were also fitted to power laws, r = cL^, using data from all 
lattice sizes. Only for Tf, the L = 6 data were discarded. The fit results are 
summarized in the following table: 





c 


z 


xVdof 




1.04(2) 


0.035(7) 


0.90 




1.30(3) 


-0.044(7) 


0.57 


n 


0.60(2) 


-0.028(8) 


0.71 



These numbers should be compared with the corresponding ones from the 
single cluster algorithm [|15]: 



L 


Te 


T^ 


16 


1.36(2) 


1.01(2) 


24 


1.50(3) 


1.06(2) 


32 


1.72(4) 


1.14(3) 


48 


1.90(6) 


1.20(4) 


64 


1.97(5) 


1.20(3) 



These r's are fitted by te oc and oc We conclude that 

the exponents z of the wall algorithm are smaller than those of the single 
cluster algorithm. For the lattice sizes in question, however, the actual r's 
are of similar size. 

For the measurements of Za/Zp we employed a variant of the boundary 
flip algorithm [^], where the cluster that has been built is not flipped but 
used to construct an observable. 

For details on the statistics and lattice sizes of all the runs, see table ^. 
The total CPU consumption of the standard action runs was about 1.1 years 
on a 200 MHz Pentium Pro PC for the wall cluster simulations. All simu- 
lations with the demon algorithm for the standard Ising model took about 
half a year on a 200 MHz Pentium Pro PC. 

3.2 Improved Action 

A first estimate for jSc was obtained by locating the crossings of Ri{L) 
with i?j(2L). We used lattices of size L = 4,8,16,32 and obtained Pc = 
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Standard Action 


Improved Action 




Cluster 






Multispin 




Cluster 


L 


stat /3000 


k 


L 


stat/32 • 10^ 


L 


stat/ 10^ 


4 


100000 


4 


2 


200 


4 


30 


6 


100000 


4 


3 


200 


6 


60 


8 


100000 

A \J \J \J \J \J 


4 


4 


200 


8 


36 


10 


50100 


4 


5 


200 


10 


10 


12 


33000 


7 


6 


200 


12 


10 


14 


25911 


4 


7 


195 


14 


10 


16 


18207 


7 


8 


155 


16 


13 


20 


10950 


4 


9 


120 


18 


8 


24 


9094 


7 


10 


105 


20 


3 


28 


7504 


7 


11 


105 


22 


7 


32 


4811 


7 


12 


100 


24 


3 


40 


3105 


7 


13 


56 


28 


6 


48 


2620 


7 


14 


51 


32 


6 


56 


1501 


7 


15 


50 


36 


5 


64 


1536 


7 


16 


50 


40 


7 


80 


703 


7 


17 


50 


48 


5 


96 


506 


10 


18 


44 


56 


2 


112 


205 


10 


19 


44 






128 


203 


10 











Table 3: Left part of table: Statistics of the runs with stan- 
dard action aX (3 — 0.2216545. stat denotes the number of mea- 
surements. Two measurements were always separated by k up- 
dates with the wall cluster algorithm. Middle: Statistics of the 
multispin-coding runs. Per measurement the demons are once 
refreshed canonically, and there are 8 microcanonical updates of 
the spin-system. Last two columns: Number of measurements at 
D) = (0.383245, 0.624235) in case of the spin-1 model. Two 
measurements were separated by three cluster updates and one 
Metropolis sweep. 
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0.383245(10), with D given by eq. i.e. = 0.624235. Note that this 

was still a preliminary estimate to be refined later. 

Monte Carlo simulations were then performed at jS = 0.383245, fixing D 
according to eq. p^). We simulated on cubic lattices with linear extension L 
ranging from 4 to 56, using the single cluster algorithm in alternation with 
a Metropolis procedure to maintain ergodicity. Two measurements were 
separated by three single cluster updates and one Metropolis sweep. As for 
the standard action we used clusters built with the boundary fiip algorithm to 
obtain estimates for Za/Zp. After each growth of the corresponding cluster, 
the work done was also exploited to perform one wall cluster step as described 
for the standard action above. More detailed information on the lattice sizes 
and the statistics can be obtained from table ^. 

The runs for the determination of (3c took about three months of CPU on 
Pentium 166-MMX PCs, while the final production runs consumed a total 
of one year CPU on the same type of PC. 



4 Measured Quantities and Basic Data 

The estimates for the basic quantities relevant for the present study are 
quoted in a number of tables in the Appendix. Here is an overview: 



quantities model table 



Ri, R2, Xi Enn 


Standard Action 


25| 


Ri, R2, X 


Improved Action 


26 


Enn, Ed, E 


Improved Action 


27 


dRi/dp, dR2/d(3 


Standard Action 


28| 


dRi/dp, dRi/d(3, dR2/d(3, dR2/d(3 


Improved Action 


29 



In both models we measured Ri, R2, which were already defined in eqs. 
and (H). We also measured the quantities needed to estimate the derivatives 
of Ri and R2 with respect to the couplings. Note that in case of the spin-1 
model, eq. ( ^if ) implies that 

dRi dRi dRi rri^\ 
'd^^W^ dD- ^^^^ 

For the sake of table presentation, the derivatives of the Ri were multiplied 
by a factor f{L) = /,-i/o-63^ which compensates for the leading divergent be- 
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haviour with increasing lattice size. The susceptibihties are defined through 

X = L^ (m^) , (28) 

with m being the normahzed magnetization, cf. eq. (^. In the tables, the 
susceptibilities are multiplied by a factor For the standard action we 

measured the energy per link, 

Enn = ^ II s,Sj . (29) 

<i,j> 

The same quantity was measured also for the spin-1 model. Here in addition 
we determined estimates for 

ed = ^j:-^i (30) 

and 

E = -SiENN - Ed) . (31) 

The definition of the energy E is such that it is proportional to ^ In Z, where 
the relation of P and D is given by eq. (P^. 

Standard reweighting techniques allow us to access the same set of ob- 
servables in some neighbourhood of the simulation point. 



5 Fitting the Data 

Our aim is to test the improvement which can be reached in the estimates 
for the phenomenological couplings Ri and the critical exponents using the 
improved spin-1 action instead of the standard action. We shall first describe 
our general approach to the data analysis. 

We shall present results obtained from fitting our data for the standard 
action and the improved spin-1 action to various finite size scaling laws. It 
will turn out that the estimates obtained from the standard action are always 
compatible with those extracted from the improved action. 

In order to get a first impression of the degree of improvement that can 
be obtained, the reader is invited to look directly at the tables of the Ap- 



pendix (table ^to table 29). One can easily verify the impressively different 



behaviour of the standard and improved data, respectively. 
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Q 




0.620 



20 



40 



60 



80 



100 



120 



Figure 1: R2 = Q as function of lattice size for both models. The 
upper data (stars) are from the standard action. The bars (flat data) 
belong to the improved action. Crosses are used to show standard 
action results corrected with the leading correction to scaling contri- 
bution. The final estimate obtained from the fit analysis is plotted by 
dotted lines. 
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Za/Zp 




0.540 



20 40 60 80 100 120 



Figure 2: Ri = Za/Zp as function of lattice size for both models. 
Standard action: crosses, improved action: bars. The final estimate 
obtained from our fit analysis is also plotted with an error interval. 
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1.01 




Figure 3: Rescaled (3 derivatives of Q and rescaled % at fixed Q as 
function of lattice size for both models (see text). 
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In figure 0, the Binder cumulants Q are plotted for tlie two actions, wliile 
in figure ^ tfie critical ratios Za/Zp are given (tables ^ and |26|) . 

In both figures, the final values obtained from our fit analysis are also 
given, together with error lines indicating the range of the estimated statis- 
tical plus the systematic error. Clearly the improved samples are much more 
stable and reach their asymptotic regime for very small lattice sizes, while 
the standard samples hardly get there on lattice sizes of order 10^. In figure ^ 
we have in addition plotted the standard action result with the contribution 
of the leading order correction to scaling subtracted, namely Q — 0.105 L"'''^^, 
see table ^. 

The improved action data outperform the corrected standard action data. 
This is due to the fact that our improvement procedure not only eliminates 
correction terms of the form , but also higher order corrections of the 
from L"""^, n integer, which are generated by the same scaling field. In 
particular corrections of the type L~^'^ should be present in the observables 
of the standard Ising model. Interestingly, the corresponding corrections 
have not been taken into account e.g. in the analysis by Blote et al. [0, in 
spite of the fact that 2 u; is smaller than other exponents taken into account 
in the fit ansatze. 

Because of their importance in extracting the critical exponents v and rj, 
in figure ^ we give also the derivatives of Q and the susceptibilities x- 
order to check for a residual L-dependence and to be able to compare the 
samples of the two models, the data have been rescaled. 

The Q-derivatives have been multiplied by a factor /,-i/o-63 ^ggg tables ^ 
and ^9|) and normalized to their value at L = 56 and L = 128, respectively, 
for the improved and standard action. 

Anticipating some results which will be given below, the x's have been 
transformed by x ^ (x ~ c) L^~^c?~^, where c, rj, and d are fit parameters. 



The values of the various amplitudes were taken from tables ^ and [19 



The authors of ref. p, |^ claim that the apparent improvement in the 
scaling behaviour that we demonstrated above does not lead to reduced er- 
rors in final results for critical exponents. Taking the data of the improved 
model by itself they are in fact right. Since the coupling parameters of the 
model are determined numerically one has to expect some small residual L""^ 
corrections, which lead to systematic errors when not taken into account in 
the fit ansatze. 

However, we shall demonstrate that it is well possible to estimate the 
effects of residual corrections to scaling in a systematic way. We shall exploit 
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the fact that ratios of correction to scahng amphtudes are universal. Given 
the parametrization 



Ri{L)=R* + r,L-^ + 



(32) 



and 



dRi 



(33) 



the ratios rj/r^ and hi/vj are universal, i.e., do not depend on the details of 



the models chosen |T^. In particular, they are the same for the standard 
and improved action. These ratios can be obtained from the analysis of the 
standard Ising model data and then used in order to estimate the residual 
correction to scaling amplitudes in the improved action results. 



5.1 i?2 at fixed Ri 

Analysing the standard model data, it turns out that the Binder cumulant 
Q evaluated at a fixed value of Za/Zp is the optimal (at least concerning 
the observables that we measured) quantity to detect corrections to scaling.[f| 
Optimal means here that the relative statistical error of the leading correction 
to scaling amplitude is the smallest. We computed Q at Za/Zp = 0.5425. 
This means that first the f3{L) is computed where Za/Zp takes the value 
0.5425. Then Q is computed for these I3{L). In the following we use Q to 
denote this quantity. In principle Za/Zp could be fixed to any value, however 
for practical reasons it is preferable to take a good approximation of {Za/Zp)* . 
To leading order Q should behave as 

Q = g* + rL~^ + . . . . (34) 

First we fitted the data obtained with the improved action. We used as 
input uj = 0.81 from ref. |Tl[]. The results are given in table ^. We see that 
there is still a small amplitude for corrections to scaling present in the data. 
The value of D where the leading order corrections to scaling exactly vanish 
should be slightly larger than the one used in this study. 

In order to quantify the improvement that is achieved we have to compare 
with the data from the standard Ising model. To obtain a consistent result 
of Q* from small lattices we had to include a subleading correction to scaling 
term r' L'"^^ in the ansatz. The results are given in table |. 

T Doing finite size scaling with quantities taken at /3- values where a phenomcnological 
coupling is kept to a fixed value is inspired by ref. |]l8| . 
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Improved Action 





Q* 


r 


xVdof 


6 


0.62369(13) 


0.0039(10) 


1.16 


8 


0.62369(12) 


0.0040(11) 


1.24 


10 


0.62362(14) 


0.0047(13) 


1.29 


12 


0.62371(16) 


0.0036(16) 


0.50 



Table 4: Fit of Q with correction to scaling term 



Standard Action 



-^min 


Q* 


r 


r' 




6 


0.62294(8) 


0.1191(11) 


0.0574(33) 


1.73 


8 


0.62326(12) 


0.1131(19) 


0.0825(67) 


0.76 


10 


0.62343(14) 


0.1092(27) 


0.102(12) 


0.60 


12 


0.62329(18) 


0.1128(43) 


0.081(24) 


0.56 



Table 5: Fit of Q with correction to scaling terms 



Starting from Lmin = 8 the fits have a small x^/dof. The result for Q* 
obtained from Lmin = 10 is consistent with the results from the improved 
action. There is a clear signal for the leading order corrections to scaling. 
The corresponding amplitude is stable when Lmin is varied. Also subleading 
corrections are well visible. 

Concerning the comparison of the two models we conclude that leading 
corrections to scaling in the improved model are reduced by a factor of about 
0.11/0.004 = 28 compared with the standard Ising model. In order to com- 
pute systematic errors due to neglecting L~'^ corrections in the analysis of 
the data obtained from the improved model we assume (taking into account 
the errors in the amplitudes) that leading corrections to scaling are reduced 
at least by a factor of 22 compared with the standard Ising model. The 
reduction from 28 to 22 takes into account that we know the ratio of the r's 
only up to some statistical error. 

Using the universality of ratios of corrections to scaling amplitudes this 
reduction has to be the same for all quantities. Hence we can take the 
correction amplitudes obtained form the standard action and divide it by 22 
to obtain a bound on the leading order corrections that are to be expected 
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Standard Action 



r . 


^1 


B 


' 1 


/dnf 

A / UUi 


10 


0.54275(14){49} 


0.22165446(16){21} 


0.0347(10) 


1.81 


14 


0.54304(23){44} 


0.22165431(19){15}* 


0.0318(22) 


1.88 


16 


0.54334(26){27}* 


0.22165417(18){11} 


0.0281(26) 


1.69 


20 


0.54269(40){24} 


0.22165443(21){7} 


0.0369(50) 


1.43 




K2 


/3c 


r2 




12 


0.62201(6){79} 


0.22165363(16){64} 


0.11166(44) 


1.30 


16 


0.62244(16){40} 


0.22165405(23){25}* 


0.1077(15) 


1.01 


20 


0.62292(31){21}* 


0.22165440(28){9} 


0.1017(38) 


1.09 


28 


0.62321(61){12} 


0.22165457(40){4} 


0.0973(87) 


1.33 



Table 6: Fitting separately the i?j with eq. (p5|), fixing oj = 0.81. 



in the case of the improved modeL 



5.2 Fitting Ri and R2 



We first fitted the Ri in order to obtain estimates for the phenomenological 
couphngs i?*, and, in addition, for the critical coupling j3c. Here and in the 
following we use as an estimate of the leading correction to scaling exponent 
u = 0.81(2) from ref. ill. 



5.2.1 Standard Action, Fit Ri 

In the case of the standard action we fitted our data with the ansatz 

Ri{L, (3uc) = R* + ^{L, Pmc) A/3 + r, L"- . (35) 

In addition to a correction to scaling term of the form L"*^ we included a 
term which (to first order) corrects for deviations from being at criticality. 
A/5 is the difference between the critical jSc and (3 = 0.2216545. Fitting Ri 
to this law, we first fixed u = 0.81. The fits were done on a sequence of 
data sets obtained by discarding data with L < Lmin- The results for the 
fit parameters as function of the smallest lattice size included are given in 
table |. The procedure used to compute the estimates of the systematic 
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Standard Action 





R*2 


r2 


UJ 






6 


0.623670(82) 


0.13732(51) 


0.9384(37) 


0.22165434(19) 


2.60 


8 


0.62433(14) 


0.1449(18) 


0.9802(93) 


0.22165489(21) 


0.89 


10 


0.62455(22) 


0.1490(40) 


0.998(18) 


0.22165503(22) 


0.94 


16 


0.62494(77) 


0.174(20) 


1.070(77) 


0.22165515(43) 


0.92 


20 


0.6252(22) 


0.22(11) 


1.16(25) 


0.22165517(61) 


1.12 



Table 7: Fit of R2 with eq. (H). 



errors (curly brackets) is discussed at the end of this section. In the table we 
mark with an asterix the value of Lmin where the systematic error estimate 
becomes equal or smaller than the statistical estimate. 

The fits are reasonably stable. Figure ^ shows the fit results for Rl and 
/3c as function of Lmin. In addition, the figure gives an impression on the 
dependence of the estimates on the choice of the fixed uj, by varying its value 
through 0.78 to 0.83. The conclusion is that the dependence on the choice 
of LV is negligible compared to the statistical errors and the systematic errors 
quoted in table ^ Note that here and in the following all error ranges for 
final fit results in the plots (usually indicated by a set of horizontal lines) are 
obtained by adding the statistical and systematic errors. 

Let us now turn to the cumulant R2. The results are given in table ^ 
and in figure ^ In the figure also the results of table where we have let 00 
a free parameter, are shown. The outcome of uj is significantly larger than 
0.81, which we presently accept as a reliable estimate. We interpret this as 
a sign that this fit parameter tries to compensate for the lack of even higher 
correction exponents. 

We also fitted the two Ri together. We did this by allowing u to vary 
freely. For Ri, the exponent u is the only correction to scaling, whereas for 
R2 we included additionally a correction term with exponent x. We checked 
the two possibilities to either fix the exponent x to some value or subject it 
to fitting in order to assume some effective value. Indeed, beyond the leading 
correction to scaling exponent u, there are several exponents that could enter 
the game, like for instance 2uj, x ~ 2, or 1/z/ + tu. In table ^ we present the 
fit results. 

Notice that the amplitude of the leading correction to scaling is much 
smaller for Ri, as expected, and for R2 also the amplitude of the next- 
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0.5445 



Za/Zp*, Standard 




Pc, Standard 



0.2216548 




Figure 4: Rl and /3c, standard action, from fits to eq. (35), with 
fixed exponents uj = 0.78, 0.81, 0.83. 
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0.626 



Q*, Standard 




0.620 



Pc, Standard 




Figure 5: R2 and /3c, standard action, from fits with eq. (35). Dia- 
monds: uj = 0.81 fixed. Bars: lo free fit parameter. 
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Standard Action 



X 


Rl 


n* 
K2 






1.62 

2.0 

2.4 

2.30(41) 


0.54120(19) 
0.54256(15) 
0.54305(13) 
0.54285(31) 


0.62295(25) 
0.62261(25) 
0.62325(24) 
0.62261(29) 


0.22165516(11) 
0.22165443(11) 
0.22165437(11) 
0.22165428(14) 


0.670(37) 
0.772(31) 
0.865(27) 
0.800(58) 


final 


0.54334(26){27} 


0.62292(31){21} 


0.22165431(19){15} 





X 


n 


r2 


r'2 


xVdof 


1.62 


0.0324(20) 


0.058(11) 


0.168(19) 


4.40 


2.0 


0.0330(19) 


0.0907(87) 


0.152(25) 


1.60 


2.4 


0.0373(17) 


0.1151(72) 


0.136(36) 


1.31 


2.30(41) 


0.0334(28) 


0.100(32) 


0.184(25) 


1.63 



Table 8: Fitting together Ri and R2 for Lj^in = 8. For Ri we 
used eq. (|35|), while for R2 also an effective next-to- leading correction 
r2 was added. 



to-leading correction to scaling ri^ is significantly different from zero. The 
effective next-to-leading exponent is of order two, and actually compatible 
with that appearing in the improved case. This could well be a coincidence. 

Finally, the errors of this fit are smaller than the (safer) errors obtained 
from the fits with fixed u. However, as a consequence of the limited number 
of lattice sizes available, the fit with many parameters can not be checked 
for stability with respect to varying L^^in- 

In order to get a handle on the systematic errors due to the neglection of 
higher order corrections, we used the result from the previous simultaneous 
fit of Rl and R2 to estimate the effective next-to-leading exponents. We 
assume this exponent x to be of order 2. Notice that from table |^ it turns 
out that the value of r'2 does not depend strongly on x. The same holds for 
r[, which is estimated to be of order 0.02 by a 4-parameter fit with u = 0.81, 
X = 2.3, on Za/Zp. 

Then we fitted separately to eq. (|35| ) the quantities defined by 



Ri{L) = R,{L) - r',L-^ , (36) 

where Ri{L) are the original (Monte Carlo) data, and r- and x have fixed 
values determined by the fits discussed above, namely r[ = 0.02, r'2 = 0.184, 
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and X = 2.3. 

Roughly speaking, the Ri{L) are the original data after subtraction of an 
estimate of the subleading correction to scaling contamination. 

The absolute values of the differences between the R* and [3c obtained in 
this way and those obtained from fitting Ri{L) with the same equation are 
the estimates of the systematic errors given (inside curly brackets) in table |^. 

In summary, from table ^ we obtain the following final estimates (labelled 
with a * in the table): 

Rl = 0.54334(26){27} 
standard action: R*^ = 0.62292(31){21} 

Pc = 0.22165431(19){15} from R^ 
Pc = 0.22165405(23){25} from R2 

The final estimates together with their error-lines are also given in fig- 
ures I and Ij. 



5.2.2 Improved Action, Fit Ri 

We fitted the data to 

R,m Puc) = R* + -^{L, Puc) AP . (37) 

Again we included a term which (to first order) corrects for deviations from 
being at criticality. Pmc is our simulation coupling 0.383245, and the dRi/dp 
are taken from table |29|. The fit parameters are R* and Pc-, entering through 
A/5 = Pmc~(^c- We first fitted separately Ri and R2 in order to compare their 
scaling behaviour. The results are reported in table ^. For both quantities 
the fit parameter estimates are quite stable. 

Here and in the following, for what concerns estimates obtained from 
the improved action, in addition to the usual error bars (quoted in usual 
brackets), we quote for each fit parameter two systematic errors. The sys- 
tematic errors are meant as an estimate of the uncertainty due to corrections 
to scaling terms. The first one, square brackets, estimates the error made 
neglecting a leading correction scaling term. The second one, curly brackets, 
as in the standard action case, estimates the error made neglecting higher 
subleading corrections to scaling. They were obtained in a well defined way 
to be described at the end of this section. 
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Improved Action 



-^min 


Rl 




xVdof 


8 


0.54213(8) [24] {66} 


0.3832470(8) [8]{36} 


1.78 


10 


0.54240(10)[19]{32} 


0.3832453(9) [5]{15} 


0.79 


12 


0.54251(11)[17]{21} 


0.3832447(9) [4] {9}* 


0.49 


16 


0.54260(15)[16]{15}* 


0.3832442(11) [4] {5} 


0.46 


20 


0.54252(25)[13]{10} 


0.3832446(14) [2] {3} 


0.55 




r?* 


(3c 


xVdof 


8 


0.62447(6) [76] {30} 


0.3832499(10) [46] {29} 


2.22 


10 


0.62429(7)[62]{15} 


0.3832479(12)[31]{12} 


1.39 


12 


0.62414(8)[57]{11} 


0.3832465(11)[26]{9} 


0.42 


16 


0.62405(12)[50]{7} 


0.3832457(14) [20] {5} 


0.31 


20 


0.62393(18)[43]{4} 


0.3832447(18) [17] {3}* 


0.28 



Table 9: Fitting separately the Ri with eq. (|37|). The numbers in 
square and curly brackets are estimates of the systematic errors (see 
text). 



Improved Action 



-^min 


Rl 


E?* 




xVdof 


8 


0.54206(7)[11]{65} 


0.62441(5) [64] {32} 


0.3832481(6)[13]{34} 


2.19 


10 


0.54231(9) [7]{32} 


0.62421(6)[51]{16} 


0.3832463(7) [8] {14} 


1.21 


12 


0.54245(11)[7]{22} 


0.62408(7)[46]{11} 


0.3832453(8) [6] {9} 


0.51 


16 


0.54254(14)[6]{14}* 


0.62399(8) [40] {7} 


0.3832448(9) [6] {5}* 


0.42 


20 


0.54252(22) [4] {9} 


0.62393(13)[35]{5} 


0.3832446(11) [4]{3} 


0.41 



Table 10: Fitting simultaneously the Ri with eq. (|37|) . 



30 



Improved Action 



X 


Rl 


K2 






0.81 


0.54423(23) 


0.62314(13) 


0.383241(1) 


1.52 


2.0 


0.54287(12) 


0.62394(7) 


0.383244(1) 


1.00 


2.5 


0.54268(11) 


0.62405(7) 


0.3832446(8) 


0.96 


2.46(56) 


0.54269(17) 


0.62405(13) 


0.3832446(11) 


1.00 


final 


0.54254(14) [6]{14} 


0.62393(13)[35]{5} 


0.3832448(9) [6] {5} 





X 


n 


r2 


0.81 


-0.014(1) 


0.0076(6) 


2.0 


-0.07(6) 


0.036(4) 


2.5 


-0.162(13) 


0.080(8) 


2.46(56) 


-0.15(57) 


0.07(22) 



Table 11: Fitting simultaneously the Ri with eq. (|3^), see text. 



We also fitted all the -Rj-data together with three parameters {Rl, R2, 
and Pc)- The results are presented in table |1^. 

Before describing the procedure we used to estimate systematic errors for 
the improved action case, let us report the results from a simultaneous fit of 
Rl and R2 with the law 

R,{L, /3mc) = R* + -^{L, Puc) A/3 + r, L"- . (38) 

As already discussed in the previous section, x represents an effective expo- 
nent. 

Our fit results are summarized in table 0. In the first three lines fits 
with fixed x are shown, while in the fourth x is a free parameter. All these 
fits correspond to L^i^ = 6. In the fifth line our final estimates are given for 
comparison. 

If we force the x-exponent to assume the leading correction value u = 
0.81, we observe that x^/dof is a little larger than for the other values of 
X. The correction amplitudes rj are very small. The main problem of this 
fit is that the ratio ri/r2 is completely inconsistent with that found for the 
standard Ising action. 

Leaving x free, it tends to choose a value around 2.5. However the corre- 
sponding fit does not improve on the estimates of the fit without correction 
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0.6250 



0.6245 



0.6240 



0.6235 



0.6230 



Q*, Improved 




Pc, Standard 




Figure 6: i?2 and /?c as function of smallest lattice size included in 
a common fit of Ri and R2 data to eq. (|37|). 
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to scaling, though consistent with it. Compare also the fits with x fixed to 2 
and to 2.5. 

Note that the x^/dof values do not allow to discriminate between the 
different solutions, i.e. the different exponents. Obtaining estimates from a 
(relatively) large lattice limit of fits without correction to scaling is, in this 
case, a safer procedure compared to a multi-parameter fit on all lattice sizes, 
in principle legitimate but in practice difficult to control. 

Being the x^/dof not a good signal of the presence of correction to scal- 
ing terms, we estimated the systematic error due to neglecting the leading 
correction to scaling term as well as that due to subleading corrections. 

For the estimates of errors due to leading corrections to scaling we used 
the following procedure (systematically used also in the following sections). 
From table |6| we know with good precision the leading correction amplitudes 

(S) 

Tj for the standard action. From the universality argument discussed in 
section 5.1, we assume that the corresponding amplitudes for the improved 
action are given by rp^/22. Therefore, analogously to eq. (|36D , we define the 
tilde quantities 

JS) 

R,{L) = R,{L) ~ ^ . (39) 



As in the previous section, fitting Ri and Ri with eq. { p7\j and taking the ab- 
solute value of the differences of the outcome parameters gives the estimates 
reported in square brackets. 



The information from the fits of table |TT] with the extra exponent x can 
instead be used to estimate in a systematic way the effect of ignoring the 
corresponding corrections. We define again 

R,{L) = R,{L)-nL-\ (40) 



where x and are parameters obtained from the fits reported in table |TT 
namely ri = 0.07, r2 = —0.15, and x = 2.46. Repeating the steps followed 
above, one obtains the estimates given in curly brackets. 

Using the fit estimates marked with an asterix, i.e., where statistical 
and systematic error estimates are of the same order, we obtain our final 
estimates: 

RI = 0.54254(14) [6] {14} 
improved action: = 0.62393(13) [35] {5} 

= 0.3832448(9) [6]{5} 
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Standard Action 



Lmin 


ai 


V 




10 


-1.4723(7) 


0.62981(7) [102] 


5.38 


20 


-1.4798(26) 


0.63045(21) [85] 


1.00 


28 


-1.4747(36) 


0.63008(27) [69] 


0.88 


40 


-1.4757(79) 


0.63014(54)[55]* 


0.96 


48 


-1.4699(13) 


0.62977(88) [48] 


1.09 






V 




14 


0.87050 (71) 


0.6331(1)[31] 


3.35 


24 


0.8592 (29) 


0.6315(4) [18] 


0.66 


40 


0.8501(76) 


0.6305(9) [13] 


0.26 


48 


0.850(12) 


0.6304(13)[11]* 


0.31 


56 


0.848(14) 


0.6303(15)[11] 


0.39 



Table 12: Fitting dRi/df5 (top) and dR2/d(3 (bottom) with eq. (|T|). 



The i?2 and Pc (as function of Lmin) are shown in figure ||. In the figures, 
our final estimates are also indicated by dashed lines, with the error intervals 
bounded by dotted lines. 

5.3 Fitting the Derivatives of the Ri 

5.3.1 Standard Action, Fit dRi/d(3 

We first fit the dRi/dp without correction to scaling, 

^ = a,LV^ (41) 

The corresponding results are summarized in table As expected, both 
quantities suffer from strong corrections to scaling. Let us first estimate the 
systematic error due to the leading correction. We followed the procedure 
described in section 5.2. Namely, we made fits with the ansatz 

^ = a,LV-(l + 6,L--), (42) 
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Standard Action 





02 


V 


62 




8 


0.8418(20) 


0.62999(25){93}(20) 


0.1103(59) 


0.48 


10 


0.8396(29) 


0.62973(36){61}(17) 


0.1181(95) 


0.47 


12 


0.8395(36) 


0.62973(43){46}(14)* 


0.119(13) 


0.52 


14 


0.8333(48) 


0.62908(55){37}(16) 


0.147(19) 


0.35 



Table 13: Fitting di^a/d/? with eq. (H) and fixed = 0.81. 



then we defined 

d(5 d(3 ' ' ' ^ ^ 

and finally we fitted the tilde quantities to eq. (|4l|) , fixing oj = 0.81. The 



differences in the u exponents are given in the square brackes of table \L2. 



The leading amplitude corrections 6j can be found in table |T3| and table |14 
The derivative of R2 suffers from stronger systematic effects than the Ri 
derivative. Therefore we include the u correction into the fit ansatz (namely, 
we fit with eq. (^2])) and compute the systematic error made neglecting fur- 
ther subleading correction to scaling. The fit results are given in table ^ and 
in figure 0, where also the zz-exponents obtained from the fit with eq. (?T) 
are reported for comparison. 

In the table we have included as a third error bar (in () brackets) esti- 
mates of the systematic effect from varying u from 0.77 through 0.85. This 
covers a 2-(j-interval around the uj value 0.81(2) that we decide to use. Again, 
the systematic error estimates in curly brackets take into account the omis- 
sion of next-to-leading corrections to scaling. They were computed with the 
procedure discussed at the end of section 5.2.1. We fitted with the ansatz 

f) n 

^ = a,L'/^{l + hL-- + b[L~^), (44) 

and defined the tilde quantities subtracting a contribution 0262-^^/^"^, with 
an estimate 62 = —0.1. 

In table ffl we give for comparison the results of fitting i?i to eq. (H3) for 



-^min = 8. Notice that in this case the data signals the presence of a stronger 
next-to-leading correction to scaling. The u obtained in this way is anyhow 
consistent with the final estimate obtained from the R2 derivative. 
We quote as our final estimates 
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Standard Action 



X 




V 


hi 




xVdof 


1.62 


-1.446(12) 


0.62873(69) 


0.175(44) 


-0.65(11) 


0.94 


2.0 


-1.460(10) 


0.62943(61) 


0.091(32) 


-0.83(13) 


1.0 


2.4 


-1.469(9) 


0.62993(55) 


0.046(26) 


-1.27(20) 


1.1 



Table 14: Fitting dRi/d(3 with eq. (||) and fixed u = 0.81 for 



0.638 



, Standard 




0.624 



Figure 7: v resulting from fitting dR2/d(3, standard action, without 
correction to scaling (bars) and with uj = 0.81 (diamonds). 
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standard action: v = 0.63014(54) [55] (from table [ID 



0.62973(43){46}(14) (from table [TD 



5.3.2 Improved Action, Fit dRi/djd 

We fitted our data for the derivatives of the Ri with respect to /3, according 
to eq. (0). The results are given in table |T5[ The i/-estimates of this table 
are plotted in figure ^j. 

Obviously, the derivatives of the cumulant scale better than those of 
Za/Zp\ while the cumulant's derivative gives a small x^/dof already for 
Lmin = 6, for the i?i's derivative one needs Lmm = 18 in order to have a 
small x^/dof and to reach stability of the result. However, also in this case 
the leading correction to scaling is strongly suppressed. Note that the range 
of lattice sizes considered here is relatively small: the difference in scaling be- 
haviour is actually due to a bigger amplitude of the next-to-leading correction 
to scaling, as discussed below. 

The systematic error estimates due to neglecting leading order correction 
to scaling are given, as usual, in square brackets. Their evaluation followed 
the procedure used in section 5.2.2. We defined the tilde quantities by 

d~Ri _ dRi bf^ ri/!.-o.8i /.rN 

(5) 



where Oj are given in table ^ and and the leading correction amplitudes 6, 
of the standard action are taken from table |13| and table 0. The absolute 
difference of the u obtained fitting to eq. (^) the tilde and the original Monte 
Carlo data are the error estimates. 

The systematic error estimates due to subleading corrections are given in 
curly brackets. To estimate them, we used a fit ansatz 

dRi 



8(3 



Li/^(1 + 6,L-^), (46) 



and then defined as usual the tilde quantities subtracting aibiL^^'^~'^ from the 
Monte Carlo data. The amplitudes hi are given in table |TB[ Comparing the 
fit without correction to scaling of the Monte Carlo and of the tilde data we 
obtained the estimates given in the curly brackets. 

Let us make some comments on the results given in table |16|. Including 
very small lattice sizes, the -Ri derivative's deviations from eq. (^) are strong 
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Improved Action 





ai 


V 




10 


-1.1419(12) 


0.62850(14)[32]{152} 


3.08 


16 


-1.1487(23) 


0.62924(25)[27]{74} 


0.90 


20 


-1.1543(39) 


0.62979(39)[26]{43} 


0.32 


24 


-1.1552(51) 


0.62988(51)[32]{51}* 


0.40 






V 


x7dof 


6 


0.66160(42) 


0.62969(11)[22]{41} 


0.79 


8 


0.66247(63) 


0.62987(14)[18]{26} 


0.57 


10 


0.6622(11) 


0.62982(22)[16]{15}* 


0.60 


12 


0.6626(12) 


0.62989(24)[24]{11} 


0.64 



Table 15: Fit of dR^jd^ (top) and 9i?2/5/9 (bottom) with 
eq. 



0.632 



^ , Improved 




Figure 8: Fit results for v from fitting dRi/dp, improved action, with 
eq. (^). The data with better scaling behaviour belong to dR^jdfi. 
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Improved Action 



X 


ai 


V 


hx 


xVdof 


0.81 


-1.20(56) 


0.63340(46) 


-0.163(12) 


0.68 


2.4 


-1.1561(11) 


0.62993(14) 


-1.309(31) 


0.34 


2.33(28) 


-1.1572(41) 


0.63003(40) 


-1.18(68) 


0.36 


X 




V 


h2 




2.0 


0.6637(12) 


0.63009(25) 


-0.058(29) 


0.63 


2.4 


0.6635(11) 


0.63005(23) 


-0.103(51) 


0.62 


2.8 


0.6633(10) 


0.63002(21) 


-0.190(93) 


0.61 



Table 16: Fitting dRi/d(5 (top) and dR2/d(3 (bottom) with eq. (H). 



Improved 




Figure 9: Fits of dRi/dp, improved action, with eq. (|4l|) (bars) and 
with eq. ( |4^ ) and x = 2.4 fixed (diamonds). 
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0.6290 ^ ' ' ' ' ' ' ' ' ^ 

6 8 10 12 14 16 18 20 22 24 

^min 



Figure 10: Fit results for ly from fitting dR2/dl3, improved action, 
with eq. (|4l|). The z^-estimate is given for the five /?- values 0.383243 
to 0.383247 in steps of 0.000001 (left to right). 
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enough to allow a fit with a free effective exponent, which turns out to be 
2.33(28), i.e. again of order 2. Notice again that enforcing x = 0.81 does 
not prevent the fit from giving an acceptable x^- We have to rely on the 
universality of the ratios of correction amplitudes, as discussed before, to 
rule out this fit. 

In figure |^ we plot the fits to dRi / d(3 using eq. (|^) for various values of 



L-tam with fixed exponent x = 2.4, together with the fits without correction 



to scaling of table |T^. This is sufficient to demonstrate the nice scaling 
behaviour of the derivative of Ri once an effective next-to-leading correction 
to scaling exponent of order 2 has been included in the fit. 

In the case of dR2/d(3 the corrections are too small to allow for a free 
effective exponent fit. Therefore we fix the x-values between 2.0 and 2.8, as 
suggested by the fits on dRi/djS. The //-results are nicely stable within these 
bounds and consistent with the values obtained without scaling corrections. 

These fits show clearly that the better scaling behaviour of the Binder 
cumulant derivative is due to a smaller amplitude of the next-to-leading 
corrections to scaling compared to that of the Za/Zp derivative. 

Finally, we checked also for the systematic dependence on the location 
of for i?2) which is giving the nicest result. We repeated the fits for the 
Q-derivative on data from five shifted /3-values ranging from 0.383243 to 
0.383247 in steps of 0.000001, covering thus two standard deviations around 
our (3c estimate. The results for u are shown in figure The effects of this 
variation is neglible compared with the errors of table |r^. 

We thus quote as our final estimate for u 

improved action: u = 0.62988(51)[32]{51} fromdRi/dp 

V = 0.62982(22)[16]{15} from 9i?2/<9/? 

The final estimate of u appears with dotted error-lines in figures ^ and |^. 
Recall that to this end the statistical and systematic errors were added up. 

5.4 Fitting the Susceptibility 
5.4.1 Standard Action, Fit x 

It turns out that the estimate of rj from fits of the magnetic susceptibility x 
taken at Pc depends quite strongly on the value of Pc- Taking the magnetic 
susceptibility at a fixed value of a phenomenological coupling removes this 



problem, as discussed in ref. [jT8[. One defines a function /5(L) by requiring 
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Standard Action 



-^'min 


c 


d 


V 


/ 


xVdof 


10 


-0.515(83) 


1.5600(44) 


0.03751(62) 


-0.585(12) 


0.96 


12 


-0.67(24) 


1.5548(84) 


0.0368(11) 


-0.568(27) 


0.98 


14 


-1.00(41) 


1.546(12) 


0.0358(15) 


-0.537(41) 


1.06 



Table 17: Fit of x at fixed Q = 1/1.604 = 0.62344, using eq. (0) 
with Lo = 0.81. 



Standard Action 



I-'min 


c 


d 


V 


/ 


xVdof 


8 


-1.08(8) 


1.543(4) 


0.0354(5) 


-0.096(12) 


1.22 


10 


-0.78(14) 


1.553(6) 


0.0366(8) 


-0.130(19) 


1.14 


12 


-1.2(6) 


1.544(16) 


0.0355(19) 


-0.095(61) 


1.15 


14 


-1.5(8) 


1.538(15) 


0.0348(19) 


-0.069(68) 


1.24 



Table 18: Fit of x at fixed Za/Zp = 0.5425, using eq. with 
LO = 0.81. 



that for any L the relation Ri{L, (3{L)) = const holds. The susceptibility 
is then computed at /3(L). We performed this analysis for the two cases of 
fixing Q = 0.6240 and fixing Za/Zp = 0.5425. Note that in principle any 
value for Q and Za/Zp that can be taken by the phenomenological couplings 
would work. However, for practical purposes it is the best to take good 
approximations of R*. 

We fitted our data to the ansatz: 

X iL,(3M) = c + dL'-^ {l + fL--) , (47) 

where c is the leading analytic part of X; and / L^'^ gives leading order 



corrections. The results of the fits are given in table T7 and in table IS. In 
both cases an acceptable x^/dof is reached at Lmm = 10. It is interesting to 
see that the correction to scaling amplitude / is considerably smaller in the 
case of fixed Za/Zp. Therefore it seems reasonable to assume that also L"^"^ 
corrections are smaller for fixed Za/Zp. We thus take 1] from fixed Za/Zp at 
Lmin = 10 as our final result. As estimate of the systematic error we quote 
the difference to the fixed Q result at Lmin = 10: 

standard action: r] = 0.0366(8){9} . 
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Improved Action 





c 


d 


V 


xVdof 


4 


-0.459(8) 


0.9496(9) 


0.0351(3) [17] 


0.42 


6 


-0.540(56) 


0.9516(18) 


0.0357(6) [14] 


0.28 


8 


-0.553(77) 


0.9519(20) 


0.0358(6) [14] 


0.30 


10 


-0.58(13) 


0.9525(29) 


0.0359(9) [12] 


0.32 


12 


-0.55(13) 


0.9519(28) 


0.0358(8) [12] 


0.33 



Table 19: Fitting x with eq. (H) at fixed Q = 0.6240. 



Improved Action 





c 


d 


V 


xVdof 


4 


-0.541(49) 


0.9545(17) 


0.03664(53) [20]* 


0.15 


6 


-0.538(61) 


0.9544(20) 


0.03662(64) [20] 


0.15 


8 


-0.532(62) 


0.9543(20) 


0.03657(60) [20] 


0.15 


10 


-0.511(90) 


0.9538(24) 


0.03644(75) [18] 


0.15 


12 


-0.53(16) 


0.9541(32) 


0.03651(94) [16] 


0.16 



Table 20: Fitting x with eq. (||) at fixed Za/Zp = 0.5425. 



We have checked that the uncertainty in the estimate of u leads to negligible 
errors in r]. We also performed fits of the magnetic susceptibility without 
a constant term in the ansatz. It is reassuring that the results for t] are 
consistent with those found above, when Lmin = 20 is taken. 

5.4.2 Improved Action, Fit x 

Also in the case of the improved action we computed the magnetic suscepti- 
bility at fixed Q = 0.6240 and at fixed Z^/Zp = 0.5425. We fitted our data 
with the ansatz 

x{L,/3iL))=c + dL'-\ (48) 
Here we have skipped the term L~'^. Our fit results are given in table |l^ and 



plotted in figure 11 



From the Ising model with the standard action we know that the ampli- 
tude of the correction is much smaller for x at fixed Za/ Zp than at fixed 
Q (see below). Therefore the comparison of both results gives a nice check 
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^ , Improved 




Figure 11: r/-estimates, fitting x-data at fixed Za/Zp to eq. (48), 
improved action. 



of the systematic errors introduced by the ommision of a L~'^ term in the fit 
ansatz. Since the corrections are smaller in the case of fixed Za/Zp, we quote 
the corresponding result as our final estimate. 

The systematic error due to the ommision of a L~'^ term in the ansatz is 
computed in the same way as in the previous sections. Namely, we defined 
the tilde quantities as 

x{L) = x{L)-dL-L'-^-'-'\ (49) 

where f'^^^ are taken from table |l^ for x at fixed Q and table |l8|for x a^f fixed 
Za/ Zp. Then we compare the results obtained fitting the x with eq. ( ^8] ) and 
the Xi both at fixed Q and fixed Za/Zp. The absolute differences of the r] 
obtained in this way are the estimates of the systematic errors. Notice that 
in this case the next-to-leading correction to scaling enters with an exponent 
of order 2 which is somehow already taken into account with the analytical 
contribution denoted by c in our fit ansatze. Therefore we only quote the 
systematic error due to the leading correction. 

We haven choosen the result of Lmin = 8 as our final estimate, since it is 
consistent with the result obtained from Lmi„ = 4 and Lmin = 6. Hence our 
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final estimate is 

improved action: i] = 0.0366(6) [2] . 
which is consistent with the estimate obtained from the standard action. 

5.5 Fitting the Energy 

As a consistency check we tried to obtain the exponent u from the singular 
behaviour of the energy. Since it turned out that the statistical errors of the 
results for u are larger than those obtained from the derivatives of the Ri we 
skipped the elaborate analysis of systematic errors that was used above. 



5.5.1 Improved Action, Fit E 

The expectation values of the energy E in the spin-1 model, defined through 
eq. (|3TD where subjected to a fit with 

E = aE + bEL^^''~'^ . (50) 



The results are quoted in table |21|. The dependence of the estimate for u on 
Lmin is shown in figure |12|, where the estimates obtained from the standard 
action are also given. 

The results are consistent with the previous estimate of u but with bigger 
errors. The final //-estimate, with error-lines, obtained from the cumulant 
derivative for the improved action is plotted for comparison. 



5.5.2 Standard Action, Fit E 

We first tried a fit without corrections, using eq. (|5I1|). The fit results for the 
nearest neighbour energy E^n of the standard action, defined in eq. (pOj), 



are given in table ^2 



Compared to the case of the improved action, the table starts with much 
larger Lmin. The corrections to scaling are much bigger than in the improved 
case. This becomes obvious from a look at figure |T^, which shows the depen- 
dence of the //-estimate on L^^^^. Therefore we take into account corrections 
to scaling with the ansatz 

Enn = aE + bEL A + _ (51) 



The results given in table for u = 0.81 show a nice agreement, however 
with large errors. 
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Improved Action 



^min 


aE 


hE 


V 


xVdof 


6 


1.23156(14) 


-0.9884(18) 


0.62893(39) 


0.58 


10 


1.23155(23) 


-0.9929(64) 


0.6297(11) 


0.63 


14 


1.23156(29) 


-0.987(11) 


0.6288(18) 


0.67 


18 


1.23157(47) 


-0.983(21) 


0.6281(32) 


0.77 



Table 21: Fit of the energy E with eq. (|o|). 



Standard Action 



-t'min 


aE 


hE 


V 


xVdof 


16 


0.3302003(32) 


0.7331(18) 


0.62902(36) 


0.73 


20 


0.3302028(47) 


0.7349(39) 


0.62937(72) 


0.75 


28 


0.3302008(65) 


0.7325(74) 


0.6290(13) 


0.81 


36 


0.330201(12) 


0.733(19) 


0.6290(30) 


0.97 



Table 22: Fit of the energy E^^jq with eq. (l50|) . 



Standard Action 





aE 


hE 


V 




xVdof 


8 


0.3302162(32) 


0.7693(35) 


0.63372(53) 


-0.155(10) 


1.11 


14 


0.3302057(59) 


0.754(12) 


0.6317(16) 


-0.098(44) 


0.87 


18 


0.3302093(93) 


0.756(30) 


0.6321(39) 


-0.09(14) 


0.85 



Table 23: Fit of the energy E'atat with eq. ( [5l| ) and a; = 0.81. 
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0.640 



0.635 



0.630 



0.625 



, Improved & Standard 



0.620 




4 6 



10 12 14 16 18 20 22 



Figure 12: Fits of the energy with eq. ( pO[ ) for both actions (without 
correction to scaling). 



6 Comparison with Other Estimates 

It is interesting to compare our estimates with other ones available from the 
literature. Note that extensive tables with many data from literature and 
experiment can be found in refs. 0, |Tl[], and |T9| . 

Some of the more recent estimates for the v and ?7 exponents, together 
with the critical coupling of the standard Ising model are collected in table 
24. Our estimates are given in the last two lines: the upper one gives to 
the estimates obtained with the standard action (SA) while the lower one 
those from the improved spin-1 action (lA). The underlined estimates are 
obtained by choosing our best estimate from the improved action and adding 
the statistical and systematic error estimates in order to obtain the overall 
uncertainty. 

The value we obtained for the Binder cumulant from the improved action 
Q = 0.62393(13) [35] {5} can be compared with the estimate Q = 0.6233(4) of 
ref. [0. It is reassuring that the high precision estimates of the recent years 
seem to be nicely consistent with each other. 
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Ref. 


Method 


V 


T] 


He 


i 




EPS 


0.6293(26) 


0.036(6) 








3D FT 


0.6304(13) 


0.0335(25) 




p 




3D FT 


0.6301(5) 


0.0355(9) 




M 




HT 


0.6310(5) 






pi 




MC 


0.6308(10) 






ill 


MC 


0.6301(8) 


0.037(3) 


0.2216544(3) [3] 


i 




MC 


0.6294(5) [5] 


0.0374(6) [6] 


0.22165456(15) [5] 


SA 


MC, R[ 


0.63014(54){55} 




0.22165431(19){15} 


SA 


MC, R'2 


0.62973(43){46}(14) 




0.22165405(23){25} 


SA 


MC, X 




0.0366(8){9} 




lA 


MC, R'2 


0.62982(22)[16]{15} 






lA 


MC, R'2 


0.6298(5) 






lA 


MC, R'l 


0.62988(51)[32]{51} 






lA 


MC, X 




0.0366(6) [2] 




lA 


MC, X 




0.0366(8) 





Table 24: Results of the present study from standard (SA) and 
improved (lA) actions are compared with other estimates: from 
e-expansion (EPS), field theory calculations in three dimensions 
(3D FT), high temperature expansions (HT) and Monte Carlo 
simulations (MC). The underlined estimates for the critical expo- 
nents are our best estimates together with error estimates which 
give the overall uncertainty, including systematic effects. 
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7 Conclusions 



By performing a detailed comparison with high precision resuhs of the stan- 
dard action, we have demonstrated that the spin-1 Ising model with suitably 
tuned coupling constants has remarkably improved finite size scaling proper- 
ties. We obtained estimates of very high precision for the critical exponents 
V and ?7 and two other universal quantities, the Binder cumulant Q and the 
ratio of partition functions Za/Zp. 

All estimates from the two different actions are consistent with each other. 
In spite of the higher statistics and the bigger lattice sizes of the standard 
action data, the estimates from the improved action are by far more precise. 
In particular, c.f. table 2^, the systematic errors are smaller for the improved 
action than for the standard action. 

The authors of refs. p|, ^ claim that an improvement of the action as 
discussed in this paper and in ref. does not allow for more precise estimates 
of universal quantities such as the critical exponents. In their argument they 
ignore the fact that ratios of correction amplitudes are universal. Once these 
ratios are computed for the standard Ising model, where the corrections are 
large, they allow for powerful bounds in the case of the improved action. 
E.g., leading order corrections to scaling of the Binder cumulant are much 
stronger than those of the derivative of the Binder cumulant. Therefore it is 
quite clear, that we can safely ignore corrections in the analysis of Q' 
obtain from the improved model. 

It would be worthwhile to use the present model in studies of physical 
quantities not discussed in this work and to check to what extent the im- 
proved scaling behaviour helps to get better estimates. 

An interesting question is whether and to what extent the improvement 
that we have achieved can be further enhanced. It seems tempting to study 
models with even more couplings in order to systematically reduce the effects 
from next-to-leading corrections to scaling. However, little is known about 
these higher order corrections, and one does not really know how many pa- 
rameters would be necessary to remove corrections of order with a; ~ 2. 
Note that at that level of improvement also the question of improved ob- 
servables comes into play. We thus believe that the improvement with two 
parameters which we performed is the optimal thing one can do in a system- 
atic way. Note that our improvement also eliminates subleading corrections 
of the typ , cf. the discussion in section 5. 

Last but not least, application of the ideas underlying the present analysis 
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to other models seems very promising. 
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Appendix: Basic Monte Carlo Results 



Standard Action 



L 







a/ 


Hj ^^ 7\r 
-*-'iV iV 


4 




6597860(58) 


1 324852(21) 


4310732('41'1 


5 




6537418(59) 


1 350645(22) 


4047035(33) 


6 


n 55n587('31') 


6492132(60) 


1 363074(23) 


3881589(27) 


7 




6458092(62) 


1 369365(24) 


3769937(24) 


8 


fl 54QlQ5f34'l 


6431935(72) 


1 372565(28) 


3690452(23) 


9 




6411374(85) 


1 374004(34) 


3631459(23) 


10 


548166('5n') 


6394995(96) 


1 374541(37) 


3586312(22) 


11 




6381157(99) 


1 374188(38) 


3550664(20) 


12 


54737Qf52'l 


636980(11) 


1 373593(41) 


3522077(19) 






n 636038(1 5) 


1 372799(58) 


n 349871 2(24) 


14 


n 546851(71) 


635224(17) 


1 371762(63) 


3479251(23) 


15 




0.634494(18) 


1.370545(66) 


0.3462842(22) 


16 


0.546212(76) 


0.633875(19) 


1.369375(69) 


0.3448912(21) 


17 




0.633329(20) 


1.368077(72) 


0.3436884(20) 


18 




0.632816(23) 


1.366696(80) 


0.3426440(20) 


19 




0.632386(24) 


1.365440(83) 


0.3417330(19) 


20 


0.54587(12) 


0.631947(89) 


1.36379(31) 


0.3409214(70) 


24 


0.54552(12) 


0.630659(82) 


1.35886(28) 


0.3384943(49) 


28 


0.54512(13) 


0.629767(93) 


1.35386(32) 


0.3368733(46) 


32 


0.54499(12) 


0.62906(14) 


1.34917(40) 


0.3357286(49) 


40 


0.54447(21) 


0.62805(15) 


1.34171(52) 


0.3342390(47) 


48 


0.54432(20) 


0.62720(17) 


1.33415(58) 


0.3333177(42) 


56 


0.54383(32) 


0.62705(23) 


1.32848(81) 


0.3327115(48) 


64 


0.54354(30) 


0.62690(23) 


1.32378(82) 


0.3322842(41) 


80 


0.54373(51) 


0.62568(34) 


1.3144(13) 


0.3317204(49) 


96 


0.54263(49) 


0.62585(36) 


1.3079(13) 


0.3313780(38) 


112 


0.54546(83) 


0.62412(64) 


1.2969(20) 


0.3311406(48) 


128 


0.54337(80) 


0.62518(62) 


1.2921(22) 


0.3309828(45) 



Tabic 25: Results for Ri = Za/Zp, R2 = Q, the susceptibility % 
divided by L^, and the energy per link, E^n, at /3 = 0.2216545. 
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Improved action 



L 




Q 




4 


0.53599(07) 


0.62408(05) 


0.87609(13) 


6 


0.54080(07) 


0.62491(06) 


0.88104(14) 


8 


0.54193(09) 


0.62460(07) 


0.87706(16) 


10 


0.54209(18) 


0.62464(13) 


0.87272(33) 


12 


0.54248(19) 


0.62416(14) 


0.86778(35) 


14 


0.54230(21) 


0.62426(15) 


0.86399(38) 


16 


0.54263(17) 


0.62409(13) 


0.86000(32) 


18 


0.54251(23) 


0.62410(17) 


0.85697(42) 


20 


0.54286(28) 


0.62379(21) 


0.85332(50) 


22 


0.54244(30) 


0.62404(23) 


0.85115(55) 


24 


0.54229(23) 


0.62402(17) 


0.84914(43) 


28 


0.54245(31) 


0.62386(23) 


0.84373(55) 


32 


0.54259(31) 


0.62390(23) 


0.84010(56) 


36 


0.54213(41) 


0.62420(30) 


0.83742(73) 


40 


0.54227(32) 


0.62403(23) 


0.83392(56) 


48 


0.54250(39) 


0.62395(29) 


0.82803(69) 


56 


0.54272(63) 


0.62382(47) 


0.82351(12) 



Table 26: Results for Ri = Za/Zp, R2 = Q, and the susceptibility x 
divided by L^, at /? = 0.383245. 
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Improved Action 



L 


Enn 


Ed 


E 


4 


0.283630(29) 


0.647597(17) 


1.091901(43) 


6 


0.248854(18) 


0.633033(10) 


1.152537(28) 


8 


0.233532(15) 


0.626493(8) 


1.178885(24) 


10 


0.225261(23) 


0.622948(12) 


1.193062(36) 


12 


0.220134(19) 


0.620744(9) 


1.201829(31) 


14 


0.216748(18) 


0.619288(9) 


1.207619(28) 


16 


0.214326(12) 


0.618243(6) 


1.211750(19) 


18 


0.212558(14) 


0.617481(7) 


1.214768(22) 


20 


0.211181(15) 


0.616886(7) 


1.217115(24) 


22 


0.210132(14) 


0.616434(7) 


1.218905(23) 


24 


0.209291(10) 


0.616073(5) 


1.220344(16) 


28 


0.207992(10) 


0.615514(5) 


1.222565(17) 


32 


0.207084(9) 


0.615120(4) 


1.224110(14) 


36 


0.206428(10) 


0.614838(5) 


1.225227(16) 


40 


0.205910(7) 


0.614613(3) 


1.226109(11) 


48 


0.205183(6) 


0.614299(3) 


1.227350(11) 


56 


0.204704(9) 


0.614094(4) 


1.228169(14) 



Table 27: Results for energy Enn, Ed and E, aX /3 = 0.383245. 
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Standard Action 



L 


f{L) dRi/dp 


f{L) dR2/dp 


4 




0.868988(29) 


5 




0.866166(30) 


6 


-1.45953(13) 


0.863364(31) 


7 




0.860978(32) 


8 


-1.47001(13) 


0.859086(37) 


9 




0.857510(44) 


10 


-1.47311(18) 


0.856291(50) 


11 




0.855176(52) 


12 


-1.47435(26) 


0.854265(56) 


13 




0.853443(79) 


14 


-1.47492(26) 


0.852892(87) 


15 




0.852278(93) 


16 


-1.47560(30) 


0.851801(97) 


17 




0.85112(10) 


18 




0.85076(12) 


19 




0.85033(12) 


20 


-1.47517(44) 


0.84968(52) 


24 


-1.47462(45) 


0.84937(51) 


28 


-1.47396(51) 


0.84789(58) 


32 


-1.47305(64) 


0.84710(72) 


40 


-1.47400(84) 


0.84641(94) 


48 


-1.4737(14) 


0.8462(10) 


56 


-1.4714(13) 


0.8457(14) 


64 


-1.4756(16) 


0.8454(15) 


80 


-1.4749(20) 


0.8479(21) 


96 


-1.4721(49) 


0.8451(24) 


112 


-1.4716(37) 


0.8451(41) 


128 


-1.4727(53) 


0.8433(41) 



Table 28: Results for derivatives of Ri = Za/Zp, R2 = Q, with re- 
spect to P, multiplied by the factor f{L) = at (3 = 0.2216545. 
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Improved Action 



L 


f{L)dRi/dl3 


f{L)dRi/dl3 


f{L)dR2/d(5 


/(L) dR2/dp 


4 


-1.09974(17) 


-.63210(11) 


.65894(15) 


-.09882(03) 


6 


-1.13595(17) 


-.64829(10) 


.66236(15) 


-.09685(03) 


8 


1.14621(20) 


-.65290(12) 


.66295(18) 


-.09614(03) 


10 


-1.15051(41) 


-.65498(24) 


.66276(37) 


-.09577(07) 


12 


-1.15293(45) 


-.65567(25) 


.66307(40) 


-.09571(07) 


14 


-1.15385(49) 


-.65661(27) 


.66310(46) 


-.09555(08) 


16 


-1.15432(44) 


-.65652(25) 


.66290(39) 


-.09554(07) 


18 


-1.15547(59) 


-.65692(33) 


.66294(55) 


-.09545(11) 


20 


-1.15603(71) 


-.65732(39) 


.66322(64) 


-.09553(11) 


22 


-1.15542(79) 


-.65714(43) 


.66310(74) 


-.09574(17) 


24 


1.1-5675(61) 


-.65790(34) 


.66419(55) 


-.09568(11) 


28 


-1.15565(80) 


-.65730(44) 


.66199(75) 


-.09520(20) 


32 


-1.15664(82) 


-.65793(46) 


.66283(75) 


-.09523(18) 


36 


-1.1568(12) 


-.65776(61) 


.6638(11) 


-.09525(49) 


40 


-1.15665(84) 


-.65785(47) 


.66326(84) 


-.09557(37) 


48 


-1.1563(11) 


-.65752(59) 


.6635(11) 


-.09513(58) 


56 


-1.1562(17) 


-.65739(97) 


.6628(19) 


-.0958(12) 



Table 29: Results for partial and total derivatives of multi- 
plied by the factor f{L) = at /? = 0.383245. 



56 



